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Abstract 

We use a one-dimensional random walk on D-dimensional hyper-spheres to deter- 
mine the critical behavior of statistical systems in hyper-spherical geometries. First, 
we demonstrate the properties of such a walk by studying the phase diagram of a 
percolation problem. We find a line of second and first order phase transitions sepa- 
rated by a tricritical point. Then, we analyze the adsorption-desorption transition for 
a polymer growing near the attractive boundary of a cylindrical cell membrane. We 
find that the fraction of adsorbed monomers on the boundary vanishes exponentially 
when the adsorption energy decreases towards its critical value. 
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The study of random walk models in statistical mechanics has led to a variety of practical 
applications [0 and illuminating analytical results. @ They provide many insights into the 
nature of critical phenomena observed in more complicated systems. For instance, in statis- 
tical field theory, the success of random walk models is best epitomized by its application to 
the triviality problem. |3], |J A fruitful area of application closer to the subject of this letter 
concerns the study of polymers using self-avoiding walks. || 

In this letter we propose a new model of random walks that substantially simplifies the 
study of statistical systems in arbitrary dimensions. In general, in systems with spherical 
symmetry we can decouple radial and angular motion. For such systems it is often found that 
many of the qualitative features of the dynamics are determined by the radial component 
alone. We show that nontrivial critical properties of physical quantities depending only on a 
radial separation are obtained with our model of random walks on concentric hyper-spheres. 

With this model we also extend the powerful techniques developed for random walk 
models to access the plethora of phenomena that spatial boundary conditions induce. It is 
well known that the introduction of a boundary can alter the physical character of a statistical 
system. [[/]] As an example of the subtle effects that our method reveals for statistical systems 
near curved boundaries, we discuss the adsorption-desorption transition of polymers growing 
near an attractive boundary. H In the limit of an infinitely extended polymer, a finite fraction 
P(k) of monomers gets adsorbed on the boundary as soon as the attractive potential k on 
the boundary increases above a critical value. For planar boundaries one generically finds 
that P{k) vanishes linearly when k approaches k c in the adsorbed phase. With our model we 
can study such a system near a curved boundary like a cylindrical cell membrane. |]] Here, 
the configurational entropy due to the "open space" available to the polymer chain is bigger 
than that in the neighborhood of a flat boundary. Accordingly, one would expect that the 
transition is weaker. In fact, with our model we show that for a cylinder of radius m > in 
monomer units and for An = k — k c — > + , the asymptotic behavior of P(k) is given by 
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P(K) ~ 81 (m+W (1) 

This asymptotic expression ceases to be valid when Ak ~ 0.188/ (m + 1) where we observe 
a crossover to linear growth in the exact expression of P(k). 

In this letter we first describe our model of random walks on a hyper-spherical geometry. 
We calculate the phase diagram for a simple percolation model to demonstrate the effects 
that random walks in such a geometry bring about. Then, we present our calculation of the 
adsorption transition for the case of a cylindrical geometry which leads to Eq. (1). 

Consider an infinite set of concentric and equally spaced spheres in arbitrary spatial 
dimension D. S n , n — 0, 1, 2, 3,..., designates the surface of the nth sphere from the 
center with area 2 r K D / 2 n D ~ l /T(D /2). We define a random walk on such a configuration in 



the following way: ||10|| Let c n< t, t > 0, n > 1, be the probability for a random walker to be 
located anywhere between S n ^i and S n at time step t. In this way, the angular position of 
the walker with respect to the origin is completely averaged away. The walker may have 
the choice to stay inside the region between S n -i and S n with probability P s tay(n), walk 
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outward with a probability proportional to the total surface area of the sphere just outward, 
Pout{ n ) — n D ~ l / N{n), and similarly inward with Pj n (n) = (n — l) ^ 1 / 'jV(n), n > 2, where 
the norm M is defined through P oui + P in + P sta2/ = 1. Thus, we find a 1 + 1-dimensional 
evolution equation to describe the behavior of the walker in this P-dimensional geometry: 

C n ,t = Pout{n ~ 1) Cn-l,t-l + Pstay{n) C n , t _i + P in (n + 1) Cn +1)t _i. (2) 

To be consistent with Ref. 10, we designate n = 1 to be the innermost region, leading 
to the boundary condition P. in {l) = (which is implicit for D > 1). In Ref. 10, the 
properties of Eq. (2) were analyzed for the initial condition c n) o = 5 n ,o and arbitrary P out 
and Pi n with P s t ay = 0. For example, for the case of concentric spheres it was found that 
the probability of ever returning to the origin for a random walker starting at the origin 
is given by ITd = 1 — 1/C(P — 1) f° r D > 2|] and unity for D < 2. This result is in 
qualitative agreement with ITd for a hyper-cubic lattice. ]TT| In fact, it has been shown that 
one-dimensional random walks on a hyper-spherical lattice have the same scaling behavior 
as P-dimensional random walks on a hyper-cubic lattice. [|I2] 

These ideas can be extend to the following directed percolation problem: For t = 0, 
there is one "wet" site on an infinitely extended line of "dry" sites with unit spacing. At 
even times t, sites take on only integer values, at odd times t, sites take on only half-integer 
values, ±1/2, ± 3/2, . . .. If at time t — 1 two neighboring sites at i and i + 1 are "dry" , then 
at time t the site at % + 1/2 is also "dry". Furthermore, we assume that two "wet" sites at i 
and 2 + 1 at time t — 1 always produce a "wet" site at i + 1/2 at time t. The latter assumption 
makes the percolation cluster, i. e. the region of all "wet" sites at any time t, compact and 
there are only two interfaces between the compact "wet" cluster in the middle and the two 
surrounding "dry" regions. Thus, we merely have to specify the following three situations: 
On the next time step the gap between the interfaces either makes a unit step outward with 
a probability P ou t, or it makes a unit step inward with a probability Pj n , or both interfaces 
shift to the right or to the left without widening the gap with a total probability P s tay 
The behavior of such a compact percolation cluster can be mapped into the one- dimensional 
random walk in Eq. (2) by defining c nit to be the probability that the gap between both 
interfaces at time t > is of width n > 0.[13|] Note that c n ^ is now defined for all n > 0, 



implying the boundary condition Pj n (0) = here. The percolation probability P is then 
defined as the probability of the cluster to persist for all times t, i. e. n > for all t. Thus, 
P is the complement of the probability of ever returning to the origin for a random walker in 
Eq. (2). Essam[TJ|] has analyzed this problem in depth for the case P out = q 2 , Pi n = (1 — q) 2 , 



and P s tay = 2g(l — q) where q, < q < 1, is a constant independent of the width of the 
cluster. Generally, we find the percolation probability [see also Eq. (3.16) in Ref. 10] 



,n=l i=l r out\ t 

reducing to Essam's result: P(q) = (2q — l)/q 2 for q > 1/2, and P(q) = for q < 1/2. 



(3) 



More interesting results can be generated from Eq. (3) by choosing distance-dependent 
coefficient functions P out and Pj n . [Note that Eq. (3) does not depend on P s t a y] For instance, 

^ refers to Riemann's £ function. 
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for P out {n) = q 2 {n + l) 5 /N{n) and P in {n) = (1 - q) 2 n s /M(n), n > 1, with P out (0) = 1, we 
obtain 5-dependent critical coefficients. We compute the percolation probability Jf5| 



oo „ 

.n=0 



n -1 



(i=«) " (n + 1) 



i^) 2 ,5,l 



(4) 



Again, for 5 = we return to Essam's result. In Fig. 1 we plot Eq. (4). For all S we 
obtain q c = 1/2. But while for 5 < 1 the transition to percolation is second order, we find 
a tricritical point at 5 = 1 and a first order transition for 5 > 1. The discontinuity is equal 
to l/((5). Such a model can be interpreted as directed compact percolation in a spherical 
configuration where n refers to the radius of the percolating bubble in dimension D — 5 + 1. 
Other choices of P ou t, Pin, and P s t ay might lead to further interesting interpretations. 

As an application of the ideas presented in the previous paragraphs, we discuss the adsorp- 
tion fraction for an extended polymer growing near an attractive cylindrical boundary. |16| 



It highlights the profound impact of a curved boundary on the critical behavior of physical 
quantities that a planar approximation could not reveal. 

In outlining the theory we extend on the beautiful treatment in Refs. 2 and 17. Consider 
a lattice consisting of an infinite set of concentric cylinders of unit spacing. Let the innermost 
cylinder - the surface of the boundary - be of integer radius m > 0, the next innermost of 
radius m + 1, and so on. Each cylinder is labeled by its radius. Now consider a random 
walk in unit steps either parallel to the length or perpendicular to these cylinders, starting 
on the boundary. A parallel step is taken with a relative weight of P s tay = 1, steps outward 
and inward are taken with relative weig hts of P out {n) = 2(n + l) s /[n s + (n + l) s ] and 
Pin( n ) — 2n 5 /[n 5 + {n + l) 5 ], respectively, for 5 = 1 and n > m,Q while on the boundary 
P out (m) = 1, Pi n (m) =0. At each step the walker picks up a statistical weight z, while 
for each step on the boundary the walker also acquires an additional weight of k > 1. For 
simplicity, we neglect self-interaction and excluded-volume effects. 

A walk with L > parallel steps has reached L levels, {/ij}^ 1 , hi > m, above or on 
the boundary. We want to restrict these walks such that \hi+i — hi\ < 1, < i < The 
transfer matrix Th i+1) hi that describes the transition of the walker from the ith to the (i + l)st 
level is given by 

T jti = Z lj ~ ll K 5 ^ [PstaySj^i + P<rat{i)Sj-l,i + P in {i)S j+X ^ ■ (5) 

The partition function Zl for walks extending L parallel steps is given by 

Z L = z L £ 6 mM T hlM T h2M . . . T hLthL _ 1 = z L b^T L e. (6) 

hi 

The total partition function for walks of all length, Z = J2T=i then evaluates to 

Z(z, k) = zb {t) T{l - zT)~ l e. (7) 

2 For 6 = 2 we can study the case of a spherical boundary while for 5 = or m — > oo we recover the case 
of a planar boundary. Negative S would correspond to a polymer growing in a cavity. 

3 This restriction on the directed walk will have no impact on the scaling behavior for L — > oo. 
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If ^max is the largest eigenvalue of T, then Z diverges for z /* z^k) = 1/A max . 

The average length of a walk and the average number of steps taken on the surface of 
the boundary are usually defined to be 

< N(z, k) >= zd z lnZ(z, k), < N s (z, k) >= k8 k In Z(z, k), (8) 

respectively. Both, < N > and < N s >, diverge for z S z^k), defining an infinitely long 
walk. We want to interpret such walks with the infinite chain limit of polymers stretched 
out along a cylindrical boundary. Then, < N s > refers to the number of monomers which 
are adsorbed on the boundary as a function of the attractive potential k. The fraction of 
adsorbed monomers P(k) is given by 

hm <N '^ K } > =- " (9) 
z/ Zoo (k) <N(z,k)> -Zoo(k) an 

Thus, 2oo( K ) marks a line in the (k, z)-plane for which P(k) is defined. To obtain a non- 
vanishing adsorption fraction it is necessary that the attractive potential k is larger than 
some critical value, k c . In terms of the eigenvalues A = X(z, k) of the transfer matrix T, k c 
is found to be the smallest value of k for which T has a bound state on the line z = z^k). 

The spectrum A of the transfer matrix T is determined by the eigenvalue problem 

{9n + z^ l g n -i + z^±|)# n+ i, n>m + 2; 

g m+1 + zg m + z 2 ^2Xf 9m+2, n = m + 1; (10) 
2(m+l) 

K,g m + KZ 2m+3 9m+l, U = 171. 

This system of equations is most conveniently examined with generating function techniques: 
Defining H(x) = Y^= m +i xn gn/(^ n + 1)) we obtain a first order linear, inhomogeneous 
differential equation for H. Bound states of T are determined by requiring both, that 
g n — > for n — > oo and that g n fulfills the boundary condition at n = m. This eigenvalue 
condition is equivalent to requiring that H(x) may have no singularity for \x\ < 1. Writing 
e = 2z/(A - 1) and 7 = (1 - v 7 ! 3 ^)^, we find 
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m+l F(|,m + l;m+|;7 2 ) 
m + 1/2 F(|,m;m+i; 7 2 ) ' 1 j 



where F is a hypergeometric function. We obtain z^k) implicitly from Eq. (11) by replacing 
A = 1/z^k). If A > l + 2z, or equivalently k > k*(z) = (l + 2z)/(l + z), there is one solution 
of Eq. (11) for A = \ max such that both conditions on the eigenvectors g are met. We obtain 
n c = 4/3, z c = 1/2 at the intersection of k*(z) and z^k). Inserting z^k) = z c — Az(k) 
and k = k c + An into the equation for z^, we get in the limit Az — > + , Ak — > + : 

1 8 

Az(k) ~ — e 9(m+i)A« > A/c< l/(m + l), (12) 
48 

neglecting exponentially smaller corrections. Using Eq. (9), we arrive at Eq. (1). In Fig. 2 
we plot the exact values of P(k) for m — 0, 1, 2, 3 and m = 00 for k > k c = 4/3. Note 
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the substantial error in estimating k c for small m using an extrapolation from values where 
P(k) appears to be linear. A measure of the overestimate on k c is given by the inflection 
point of P(k) in Eq. (1) at An = 0.188/(m + 1), beyond which the approximation in Eq. (1) 
ceases to be valid and P(k) grows linearly 

The simplicity of the dynamics in this new random walk model, and the nontrivial scaling 
obtained from it, raise interesting questions regrading the universal properties of this lattice. 
At the critical transition only a few fundamental properties of the system determine its 
behavior. In this model, the critical behavior arises from the balance between a short-range 
attractive potential and the spatial entropy. We argue that these features are sufficiently 
well represented by a random walk on a hyper-spherical lattice. It has been shown that such 



a lattice reproduces all the universal scaling properties expected of lattices. [|12j The critical 
behavior obtained on this lattice for rotationally symmetric systems should therefore reflect 
the universal critical behavior of the system. The advantage of random walks on hyper- 
spheres is to describe the critical behavior in a minimal and tractable way in comparison, 
for example, to a far more structured hyper-cubic lattice. 
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FIGURE CAPTIONS 



FIGURE 1: The percolation probability P(q, 5) for directed compact percolation in a curved 
geometry. For any S : the percolation threshold occurs at q c = 1/2, but for 5 < 1 the 
transition is second order, while for 5 > 1 the transition is first order with a discontinuity of 
P(q c ,S > 1) = l/((5), indicated by a darkened wedge. 

FIGURE 2: The exact adsorption fraction P(k) plotted for m — 0, 1, 2, 3 and m = oo and 
k > k c — 4/3. For finite m, P(k) vanishes exponentially for k — k c — > + while it crosses 
over to a linear increase when k — k c > 0.188/(m + 1). For m = oo we recover the linear 
scaling for P(k) that was found in Ref. 17. 
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